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ABSTRACT 

Accurate line lists for three molecules, BeH, MgH and CaH, in their ground elec- 
tronic states are presented. These line lists are suitable for temperatures relevant to 
exoplanetary atmospheres and cool stars (up to 2000K). A combination of empirical 
and ab initio methods is used. The rovibrational energy levels of BeH, MgH and CaH 
are computed using the programs Level and DPotFit in conjunction with 'spectro- 
scopic' potential energy curves (PECs). The PEC of BeH is taken from the literature, 
while the PECs of CaH and MgH are generated by fitting to the experimental tran- 
sition energy levels. Both spin-rotation interactions (except for BeH, for which it is 
negligible) and non-adiabatic corrections are explicitly taken into account. Accurate 
line intensities are generated using newly computed ab initio dipole moment curves 
for each molecule using high levels of theory. Full line lists of rotation- vibration transi- 
tions for 9 BeH, 24 MgH, 25 MgH, 26 MgH and 40 CaH are made available in an electronic 
form as supplementary data to this article and at www.exomol.com 
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1 INTRODUCTION 

The ExoMol project aims to provide line lists of spectroscopic transitions for key molecular species which are likely to be 
important in the atmospheres of extrasolar planets and cool stars. The aims, scope and methodology of the project are 



summarised in the previous article (Tennyson & Yurchenko 20121. In the present paper, the first in a planned series, we 



present results for the electronic ground state rotation-vibration spectra of three metal hydrides BeH, MgH and CaH, denoted 
generically as XH. The line lists are explicitly designed to extend to high temperatures. All three of these hydrides are of 
astronomical interest. 

The alkaline earth metal species BeH, MgH, and CaH share many similar properties. Early spectroscopic data on these 
three species is captured by |Huber k. Herzberg| ( [T979[ ). A summary of more recent studies, both spectroscopic and astrophysical, 
is presented for each molecule below. 



1.1 BeH 



BeH is one of simplest heteronuclear diatomic molecules, and hence it is a strong contender for being observed in contexts 
such as exoplanetary atmospheres, cool stars and in the interstellar medium. However there are only very few astrophysical 
records of BeH; for example, a detection of A 2 S 



X E emission lines of BeH in the sunspot umbra spectra by 



Wohl 



(19711 and Shanmugavel et al. (20081 



BeH has been studied extensively by quantum chemists using different ab initio methods. It is the simplest open-shell 
molecule, which made it very popular as a benchmark object for testing new ab initio approaches. Ab initio studies on BeH are 



mostly limited by the basic spectroscopic and structural constants. Pitarch-Ruiz et al. (20081 reported very comprehensive 



electronic-structure calculations on BeH with the potential energy and dipole moments curves computed for a number of 



electronic states. A very accurate ab initio ground state PEC was recently generated by Koput (20111. Le Roy et al. (2006) 



generated a 'spectroscopic' potential energy curve (PEC) of BeH obtained by fitting to the experimental values of Colin et al 
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(19831; Focsa et al. (19981; Shayesteh et al. (20031. This is currently the most accurate PEC of BeH, and we employ it here 



in our production of the line list for BeH. 

Experimental work on BeH is not so extensive: most of experimental studies of BeH have focused on electronic transitions 
( Degreef fc Colin|1974 Colin et al.|1983 Clerbaux fc Colin|1991 Focsa et al.|1998 l. Ground electronic state rotation-vibration 



spectra were reported by Shayesteh et al. (2003) 



1.2 MgH 



MgH is thought to be a likely candidate for interstellar observation but has yet to be detected (Sakamoto et al. ( 1998a )). 
However its presence in stellar spectra is well documented through observation of the A 2 U — > X 2 E + ( Sotirovski||l971 1 and 
B' 2 E + — >X 2 E + (Wallace et al. 1999 1 transitions. These transitions have been used as an indicator for the magnesium isotope 



abundances in the atmospheres of different stars (Lambert et al.|1971 Tomkin fc Lambert|1980" Lambert fc Mcwilliam] 1986 



Wallace et al. 1999 Gay & Lambert 2000 Yong et al. 2003 ) from giants to dwarfs including the Sun. The same lines of MgH were 
used to measure the temperature of stars (for example, by | Wy llcr (1961)), surface gravity (for example, by Bell et ah] ( |1985[ ); 
Bonnell & Bell (19931), stars' metal abundance (for example, by Cottrell (19781; Tomkin & Lambert (19801), gravitational 



1990), as well as for a deuterium test (Pavlenko et al 



2008) 



acceleration (for example, by Berdyugina & Savanov ( 1992 1) , to measure the temperature and pressure (Kulaczewski et al. 

MgH is an important part of stellar atmospheric models, such 



as, for example, PHOENIX ( jAllard et al.|2001[ ) which incorporates the MgH line list of |Kurucz| ( |2011 1 



MgH has been studied in laboratories since 1929 (see (Pearsc 1929)), with measurements of rovibrational transitions in 



the X E + state given by 



Zink et al. 



experimentally measured spectral lines exist for the A II 



(19901; Ziurys et al. (19931 and most extensively by Shayesteh et al. (2004a I. Many 

and B' 2 E 4 



X 2 E 4 



infrared up to v = 4. Experimental ground X 2 E + state spectra were reported by Leopold et al. ( 1986 1; 



X 2 E+ transitions ( |Bernath et al.||l985 



Shayesteh et al. 2007 Shayesteh & Bernath 20111, although the ground state transitions have only been observed in the 



Lcmoine et al. 



(1988) 



Zink et al. (19901; Ziurys et al. (19931, which also include high-resolution millimeter- wave studies of the pure rotational 



transitions with the hyperfine structure resolved. 

On the ab initio side, the calculations of the potential energy curves and dipole moment curves (DMCs) of MgH by 



Saxon et al. ( 1978 ) made a very important contribution to the development of the theory for this molecule. The results of 
this study have been used in many spectroscopic applications including the 24 MgH line list from the UGAMOP group (Week 



et al.||2003b|a I. Very recently, accurate ab initio calculations have been performed using significantly higher levels of theory 
( Mestdagh et al.|2009 Guitou et al.|2010 1, with PECs and DMCs of MgH reported. Important theoretical work for the present 
study is by Shayesteh et al. (20071, where an accurate 'spectroscopic' PEC of MgH was obtained by fitting to all existing 



experimental transitions for this molecule. We use this PEC as an initial approximation in our fits. 

The first complete MgH line list was computed by Kurucz (see |Kurucz| ( |2011p ). This line list, however, overestimates the 
opacity of MgH when included into the atmospheric models due to presence of nonexistent levels (see discussion by Week 



et al. ( 2003b l). Kurucz's line list does not include the ground rotation- vibration transitions and thus cannot be compared to 



the results of the present work. The UGAMOP line lists ( |Weck et al.||2003a|b|c \ mentioned above contain 23 315 transition 
energies and oscillator strengths over the wavelength range - 32,130 cm -1 , for all possible allowed transitions from the 
ground electronic state vibrational levels with v" ^ 11. 



1.3 CaH 



The main astrophysical interest in CaH is because of the spectra of T Tauri stars (IDi Rosa 2004). The transition lines of 



CaH were used in identifications and studies of different dwarfs ( Barbuy et al.||1993| |Reid et al.|1997 Burgasser et al 
McGovern et al.||2004 Reiners et al.||2007 l, galactic disks ( Reid et al.||1997 l, molecular clouds ( Sakamoto et al.||1998b 



2003 



and 



sun spots ( Berdyugina et al.|2006 l. One of the most important spectroscopic features of the CaH spectra of T and L dwarfs 



is the 6750 - 7050 nm band covering the A — X (0, 0) rotational tran sitions (jBurg asser e t al.|2003 Reiners et al.|[2007 k 

CaH has been studied in the laboratory since Hulthen ( 19271 and Watson & Weber ( 1935 1. Much experimental work has 



been done on the electronic transitions of CaH ( |Weck et al. 2003d| [Ram et aL 2011 Leininger & Jcung 1995] [Chen et af 



20061, but only a few studies exist on the ground state rovibrational spectrum. The most important experimental studies of 
CaH for the current work include the Fourier transform emission spectra of the E 2 Tl — X 2 E + bands ( Ram et al.|201 1 1 and of 



four ro- vibrational bands in the X ^E" 1 " state of CaH and CaD ( |Shayesteh et al.]|2004b |. These works provided very accurate 
experimental measurements of rovibronic frequencies with a very thorough analysis. Other experimental work include spectra 
and analysis of the A - X and B — X systems of CaH and CaD ( |Berg fc Klynning||T9"74"| |Martin||T98"4"| |Berg et ai]|1976 
Klynning fc Martin|1981 Gasmi et al. 1 2003 Pereira et al.|2002 |; ultraviolet-absorption spectra of the C, D, K, L, G, J, and 

Kaving & Lindgrcn 19811; magnetic 



M bands of CaH ( |Kaving et al.||1974| |Kaving fc Lindgren|l974 



1976 



Bell et al. 1979 



hyperfine structure of the A — X and B — X spectra of CaH ( Steimle et al.|| 19*87 Chen et al. | |2006[ ); and the IR ( |Petitprez 



C 
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et al.|[l989 Frum fc Pickett||1993 l and millimeter- wave ( Barclay et al.|1993 Frum et al.|1993 l transitions within the ground 
electronic state. 

From a very large number of ab initio studies of CaH we select here only most recent and comprehensive works. The 



potential energy and dipole moment curves of CaH were computed using different levels of theory (Boutalib et al. (1992) 



Leininger & Jeung (19951; Holka & Urban (20061; Kerkines & Mavridis ( 2007 1). In most cases a number of equilibrium 
constants, including r c , oj c and D c , as well as the permanent dipole moments jj, e were reported for one or several low 
lying states. To our knowledge there are no experimental absolute intensity measurements of CaH. However, a few lifetime 
measurements for a number of low lying vibrational states v = 0,1,2 of the X, A, and B states were used to estimate 
the corresponding oscillator strengths or Einstein coefficients ( Klynning et al.||1982 Berg et al.||1996 Liu et aI7||2009 1 . The 
permanent electric dipole moments of CaH were obtained experimentally for a few low electronic states using the Stark-effect 
( Steimle et al.||2004| [Chen fc SteimTe||2008 1. 

As far as the theoretical ro-vibronic studies on CaH are concerned, there have been attempts to refine PECs by fitting to 
available experimental energies or frequencies (Martin 1988. Uchara 1998. Urcna et al. 2000). However, previous ro-vibronic 
models did not include a proper description of different important couplings between electronic states. Therefore only the 



ground electronic state could be reasonably well described in this approach, for example by Uehara ( 1998 1 who also included 



the Born-Oppenheimer breakdown (BOB) effect. None of these studies take account of spin-rotational splittings. Carlsund- 



Levin et al. ( 2002 1 presented an accurate model with an explicit coupling of four lowest electronic states as well as the 



BOB effect using previously published ab initio PECs by Martin ( 1988| ). However since these PECs are of ab initio (not 
experimental) accuracy, the resulting energies were only good enough to validate the models used, but not to predict the 
corresponding frequencies of CaH and CaD. 



The most comprehensive line list for CaH was reported by 



Week et al. 



(2003(1), who produced a complete list of CaH 
970) in the frequency range 



rovibrational levels and oscillator strength covering all possible allowed and bound transitions 
up to 35000 cm -1 involving the A 2 E + , A 2 U, B/B' 2 T, + , C 2 E+, D 2 E+, and E 2 U states. This line list should be good up to 
at least 2000 K and is also included into the UGAMOP data base. In the production of this data the authors used the most 
recent set of theoretical PECs and DMCs, with adjustments to account for experimental dissociation energies and asymptotic 
limits. However the model used to compute the frequencies was not very accurate: (i) no fine structure, such as spin-rotation, 
or non-adiabatic effects, were taken into account; (ii) the potential energy curves were not adjusted to the experiment thus 
the overall accuracy is no better than 10 to 30 cm -1 . 

In this work we present detailed line lists for BeH, 24 MgH, 25 MgH, 26 MgH and 40 CaH considering only the rotation- 
vibration transitions in the X 2 T, + state. At temperatures up to about 2000 K, any population of excited electronic states 
will be insignificant compared to the population of the ground state. Taking typical values for energies in the ground and 
first excited states of XH ( |Pitarch-Ruiz et al.|2008| [Guitou et alT]|2010| [Kerkines fc Mavridis|2007| >, the thermalised, relative 
population of the lowest rovibrational level of the first excited electronic state to that of the ground state for BeH, MgH and 
CaH can be estimated to be 2 x 10 -6 , 2 x 10 -5 and 1 x 10 -4 respectively at 2000 K. In practice the excited states are all close 
to or above the dissociation limits for the species so dissociative processes will further deplete the excited state populations. 



2 METHOD 



For the species XH we solve the Schrodinger equation allowing for Born-Oppenheimer breakdown effects as given by |Shayesteh| 
et al.| ( |2007| ) using two programs, Level 8.0 ( |Le Roy|200"7) and DPotFit 1.1 ( |Le Roy||2006 l. 

We did not generate new ab initio PECs for BeH and MgH. Full sets of potential parameters representing very accurate 



PECs for these molecules were already available from Le Roy et al. (20061 and Shayesteh et al. (20071, respectively. Ab initio 



potential energy curves for CaH were generated using multiple levels of theory, beginning with a simple Hartree-Fock (HF) 
method and progressing to complete active space self consistent field (CASSCF) and multi-reference configuration interaction 
(MRC I) (|Knowles fc Werner||1988|. The best potential, with no unphysical kinks and a shape in agreement with previous 
works dKerkines fc Mavridis|2007| Shayesteh et al. 20 04b| , used the MRCI/aug-cc-pV5Z level of theory. All electronic structure 



calculations were performed with the quantum chemistry package Molpro 2010 ( Werner et al.|2010 l. 

The DMCs for all three molecules were generated ab initio employing Molpro and different levels of theory. The best 
DMCs for each molecule (shown on Figj2j) were selected based on their smoothness and agreement with published curves 
( Saxon et aT|l978| |Pitarch-Ruiz et al.||2008[ |Leininger fc Jeung||1995| |Guitou et al.||2010| >. The DMCs were calculated using 
a finite field rather than as an expectation value as this is known to give more accurate results ( |Lodi fc "T ennyson 2010[ ). 

For BeH the final DMC was computed using the aug-cc-pwCV5Z basis set and the MRCI method on a grid of 45 geometries 
distributed between 1.3 and 10 A. Our equilibrium dipole moment is —0.224 D at r c = 1.342 A. This can be compared to 
-0.199 D ( |Pitarch-Ruiz et al.||2008[ ) for the same value of r c . 

For MgH the final DMC was calculated using RCCSD(T) ( Knowles et al.|1993 20001 in conjunction with the basis sets 
cc-pCV5Z ( Woon fc Dunning|1995 \ for Mg and aug-cc-pV5Z ( Woon fc Dunning|1993[ ) for H and all electrons correlated. The 
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Figure 1. Potential energy curves of CaH (solid line), MgH (dashed line), and BeH (dotted line) in their ground electronic states. The 
associated first excited states are also indicated. 




Figure 2. The ab initio dipole moment curves of BeH, MgH and CaH in their ground electronic states selected for the line list production. 



good performance of the RHF-CCSD methods in the case of MgH was noted by Li & Paldus (19951. At r c — 1.730 A we 
obtain an equilibrium value of the dipole moment fi c of —1.371 D. This can be compared to the CCSD(T) value from Serrano 



Saxon et al 



& Canuto ( 1998), —1.356 D. Our DMC agrees surprisingly well with dipole curves computed at the SDTCI level of theory by 



(1978). 



The DMC for CaH was also calculated using RCCSD(T)/cc-pCV5Z with all electrons correlated. |Holka fc Urban| ( |2006 ) 
demonstrated importance of the electron correlation contribution (0.44 D or 17.4% of the total value), while the relativistic 
and complete basis set corrections are relatively small (about 0.009 D and 0.007 D, respectively). Our equilibrium value of the 
RCCSD(T)/cc-pCV5Z dipole moment is 2.55 D, which can be compared to the experimentally derived value of the permanent 



dipole moment 2.53(3) D measured for the ground X E 



0) state with the optical Stark effect ( Chen fc Steimle[2008 ) 



Holka & Urban (20061 recently obtained an ab initio value of 2.463 D for ^t c using a high level of theory. 



In all our calculations the ab initio DMC grid points were used directly in Level. 



2.1 Fitting the potential 



In the case of BeH we employed the spectroscopic PEC of Le Roy et al. ( 2006 1 consisting of an extended Morse oscillator 
(EMO) function corrected by BOB. We could not precisely reproduce the ro- vibrational energies reported in this paper using 
these potential parameters and Level 8.0, see columns II and III of Table [T] The deviations in the vibrational term values are 
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rather small, except for v — 11 which is about 10 cm \ therefore we decided to employ the PEC from Le Roy et al 
in the line list calculations without any additional refinements. 



(2006) 



In the case of the PEC of MgH we experienced a similar problem. Using the potential parameters from [Shayeste h et al.| 
(20071 in conjunction with the Morse Long Range potential (MLR) ( Le Roy fc Henderson|2 007) in the DPotFit calculations 
we could not reproduce frequencies obtained by Shayesteh et al. (20071 using the same program. Considering importance of 



MgH for astrophysical applications we decided to refine Shayesteh et al. (2007l's parameters through fits to the experimental 



data and using the same program (i.e. DPotFit 1.1) to guarantee the quality of the line list produced. In these fits, to avoid 
dealing with the high electronic states, we constructed a set of 'fake' transitions within the ground electronic states from 
the experimental transitions A 2 Ii — > X 2 E + and B 2 T, + — > X 2 T, + using the method of the combination differences. As 
a result our 'experimental' set included 1302 transition frequencies, 380 of which were actual experimental, while the rest 
were 'experimentally' derived values. Some of the lines, especially from highly excited vibrational excitations, had a greater 
uncertainty and therefore were given a lower weighting. Despite their relatively low quality, these data are important to fix 
the behaviour of the potential at large inter-molecular distances and hence to accurately predict the higher-energy states. The 
fit was found to be unstable unless the dissociation energy D e was held fixed - the best available value of 11104.7 cm" 1 from 



Shayesteh et al. (20071 was chosen 



The fit of the MgH potential achieved a root-mean-square (rms) error of 0.0012 cm -1 for the 380 experimental X 2 £ + — 
X 2 E + lines with a maximum difference of 0.0077 cm -1 , showing that these lines are reproduced to high accuracy despite the 
inclusion of extra lines in the fit. Compared with the experimental line list (6076 lines) calculated by combination differences, 
the rms error is 0.059 cm -1 and the maximum difference is 0.35 cm -1 . The 'experimentally' derived lines show an inherent 
uncertainty of up to about 0.2 cm -1 . 

To our knowledge there is no potential for CaH in the literature of a satisfactory accuracy. Therefore we built a new 
'spectroscopic' PEC for this molecule using the MLR potential. Fitting to the experimental data requires a good initial 
estimate for the PEC, for which our MRCI-Q/aug-cc-pV5Z potential was used. The corresponding MLR potential parameters 
were obtained through a fitting to the ab initio energies (see above) distributed over the geometry range from 0.7 to 10.0 A. 
The potential was then refined by fitting to a list of 426 transition frequencies of 40 CaH by Shayesteh et al. ( 2004b I and Ram 



et al. (2011 1. We obtained higher accuracy fits if D e was allowed to vary. Our value of D e — 15 923 cm -1 ended up somewhat 
higher than the experimental value 14 800 cm -1 ( Huber fc Herzberg||l979" l. The CaH fit achieved an rms of 0.001 cm -1 for 
426 lines, the most complete data set we could build for fitting against. It should also be noted that no differences exceeded 
0.0087 cm" 1 . 

For CaH, the most important change from the initial ab initio parameters is attributed to the equilibrium geometry. The 
initial parameters of MgH, being already empirical, however, were significantly less affected by the fitting process. 

To illustrate the quality of these PECs, Tables [I] - [3] compare our vibrational energies, T v , and rotational constants, B v , 
with the corresponding 'experimental' (when available) and other theoretical values. In case of MgH and CaH we list the 
UGAMOP theoretical values ( Week et al.|20 03a b ; Skory et al. 2003 . Wee k et al.|2003c|d I, while for BeH ab initio spectroscopic 
constants of Le Roy et al. ( 2006 I are given. When comparing to our energies it should remembered that the UGAMOP energies 



do not include the spin-rotational splitting. Our results exhibit a significant improvement relative to the UGAMOP values 
for MgH and CaH. 



3 LINE LIST CALCULATIONS 

The procedure described above was used to produce line lists for the XH species, i.e. catalogues of transition frequencies 
i>ij and Einstein coefficients Aij. The computed line lists for BeH, 24 MgH, 25 MgH, 26 MgH, and 40 CaH are given in the 
supplementary materials. The same 'spectroscopic' PEC of MgH corrected for BOB was used to produce the line lists for 
all three isotopologues of this species. The coverage in terms of the frequency range and quantum numbers is illustrated in 
Tabled 

The line list file, see TableJZ] contains transition frequencies Vij (cm" 1 ), Einstein coefficients Aij (1/s), lower (i) state term 
values E'l (cm -1 ), and upper and lower state quantum numbers for each line. Quantum numbers are defined in accordance 
with the Hund's case (b) formulation: N specifies the total angular momentum excluding spin, J = N ± 1/2 specifies the total 
angular momentum including electron spin, and e/f is the rotationless parity: the sign of the total parity is defined by (— 1)^ 



for a E + electronic state so states with J = N + 1/2 are e states and those with J = TV — 1/2 are / states (Brown et al 



19751; and the vibrational quantum number v. Since we ignore hyperfine effects due to the interaction with the nuclear spins, 
J is assumed to be a conserved quantity. The e/f labels are omitted from the BeH line list as we consider the corresponding 
energies to be degenerate. The energy levels are given relative to the corresponding zero point energy -Ejv= „ =0 . 

The spin-rotational coupling plays an important role in the case of ro-vibrational spectra both of MgH and CaH, in 
contrast to BeH, where the fine structure can be considered negligible. The program Level is not capable of treating the 
spin-rotation coupling, and therefore the corresponding ro-vibrational energies of MgH and CaH were computed with the 
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Table 1. A comparison of the theoretical and experimental spectroscopic constants of BeH in cm 



Obs. a 



Calc. 6 



Obs.-Calc. 



Le Roy et al. 1 2006 I 



Tz 

T 3 

T 4 

T 5 

T, 

T 7 

T 8 

T 3 

T 10 

Tn 

B 

Bi 

B 2 

B 3 

Bi 

Br, 

B 6 
Bi 
Ms 

Bg 

Bio 

Bu 



1986.4442 
3896.8785 
5729.2861 
7480.4528 
9145.2834 
10716.2804 
12182.323 
13525.889 
14717.860 
15709.040 

10.1648880 
9.8554335 
9.5417271 
9.220768 
8.886639 
8.533209 
8.152609 
7.729457 
7.2309 
6.5921 
5.6984 



1986.4176 
3896.8714 
5729.2618 
7480.3400 
9145.2533 
10716.6292 
12183.0505 
13526.6563 
14718.8262 
15710.4672 
16402.0939 
10.165644 
9.855734 
9.541630 
9.220113 
8.886566 
8.534221 
8.152146 
7.722324 
7.214164 
6.570639 
5.662832 



0.0266 
0.0071 
0.0243 
0.1128 
0.0301 
-0.3488 
-0.7275 
-0.7673 
-0.9662 
-1.4272 

-7.56x 
-3.01x 
0.98x 
6.55x 
0.73x 
-10.12X 
4.63x 
71.33X 
166.86X 
214.91X 
355. 38x 



10- 4 
10- 4 
io~ 4 
io- 4 
io- 4 
io- 4 
io- 4 

10" 4 
10" 4 
10~ 4 

io- 4 



1986.4163 
3896.8712 
5729.2602 
7480.3383 
9145.1320 
10716.1629 
12182.2070 
13525.7882 
14718.0819 
15709.3839 
16412.7734 
10.165715 
9.855827 
9.541718 
9.220280 
8.886738 
8.534429 
8.153772 
7.728265 
7.225306 
6.579663 
5.667413 
4.268157 



Experimentally derived by |Focsa et al.|jl998| ; b Calculated, this work. 



program DPotFit. Only the Einstein coefficients were generated using Level, i.e. without spin-rotation effects included. We 
assume that the e—f intensities are negligible, while the transition moments for e-e and /-/ transitions are equal (see also 

(1989)). This can be readily verified for MgH: the known experimental e—f transitions are 



discussion in the book by Herzberg 



indeed very weak (Shayesteh et al. 2004a I, and for CaH no such transitions were reported experimentally, probably for the 
same reason. It should be noted that the Level program computes the line strengths as 'singlet-singlet' transitions and the 
following correction factors had to be applied (Watson 2008) in order to account for the fine structure in the MgH and CaH 
transitions: 

s R (f^f) = s£ cvcl ^— } . m 

S P (f^f) = Sl evel ^—^, (2) 



2J" + 


1 — 1 


2J" + 


2' 


2.1" - 


1 


2.1" 




2J" + 


3 


2.1" + 


2' 


2.1" + 


1 



S R (e^e) = SE^Z^, (3) 
S P (e^e) = sL vel ^j^. (4) 

Here S R ^ P is the line strength of the R/P transition. 

The integrated absorption coefficient I(cJij) for a i — > j transition can be calculated for a specific temperature T using 



Z{T)u 



where gj is the degeneracy of the final state, Z(T) is the partition function, fee is the Boltzmann constant, E' is the energy of 
the higher state, E" is the energy of the lower state, u) (cm -1 ) is the transition wavenumber, and a is the factor 1.3271005 x 
10 -12 s cm -2 molecule -1 so that I has units of cm/molecule. For the X 2 E state of BeH, where the states i are doubly 
degenerate, the partition function, Z(T), is given by: 

Z(T) = g*e~ Ez/kBT = g ^ 2S + V( 2Ni + 1 )e" Bl/fcBT , (6) 

i i 

where S — 1/2 is the total electronic spin angular momentum, g ns — 8 is the nuclear statistical weight of BeH, and Ni is the 
rotational angular momentum of the nuclei of the ith state. Here we follow the HITRAN ( jSimeckova et al.|2 006) convention 



and consider the nuclear statistical weight of the molecule explicitly in Z(T). In the case of MgH and CaH in their ground 
electronic states the spin-rotational splitting is resolved and Z(T) is given by 
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Table 2. A comparison of the theoretical and experimental spectroscopic constants of 24 MgH in cm 1 . 







Pair- b 


Uun. — V_j aic . 




TTr 1 A A/TOP 


Ti 


1431.97786 


1431.9777 


0.0002 




1423.0023 


T 2 


2800.67807 


2800.6770 


0.0011 




2784.3944 


T 3 


4102.32975 


4102.3281 


0.0017 




4081.3781 


T 4 


5331.3892 


5331.3867 


0.0025 




5308.9675 


T S 


6479.6562 


6479.6547 


0.0015 




6459.3356 


T 6 


7534.8137 


7534.8076 


0.0061 




7521.0069 


T~ 


8477.9997 


8477.9959 


0.0038 




8477.2820 


T 8 


9279.6527 


9279.6436 


0.0091 




9299.8872 


Tg 


9892.7243 


9892.7201 


0.0042 




9929.2662 


Tio 


10249.4074 


10249.4050 


0.0020 




10150.2699 


Tu 


10352.25 


10352.2503 


0.00 




10231.8603 


Bo 


5.73650768 


5.7365078 


-0.01x10" 


-5 




Bi 


5.55528801 


5.55528804 


-0.003x10" 






B 2 


5.367512 


5.36756316 


-5.12x10" 


-5 




B 3 


5.169725 


5.16981489 


-8.99x10" 


-5 




Bi 


4.956539 


4.95667678 


-13.78x10" 


-5 




Br, 


4.71938 


4.71971001 


-33.00x10" 


-5 




B 6 


4.44431 


4.44473868 


-42.87x10" 


-5 




B 7 


4.1072 


4.10695883 


24.12x10" 


-5 




B 8 


3.65877 


3.65971997 


-95.00x10" 


-5 




B 9 


3.00781 


3.00747254 


33.75x10" 


-5 




Bio 


1.9687 


1.9684569 


24.31x10" 


-5 




B n 


0.884 


0.88689093 


-289.09x10" 


-5 





a Experimentally derived by jShayesteh et aL|2007|; b Calculated, this work. 



Z(T)=^ ffM (2Ji + l)e- B * /fcBT . (7) 

i 

Here the nuclear statistical weig hts g nB = 2 for 40 CaH, 24 MgH and 26 MgH, and g ns = 12 for 25 MgH. The nuclear spins of 
25 Mg and Be used to define these values are 5/2 and 3/2, respectively, whilst the nuclear spins of 24 Mg, 26 Mg, and 40 Ca are 
0. 

High resolution comparisons with specific regions of the experimental spectra for each molecule are given in Fig. [3] 
note that the experimental spectrum has been scaled to the theoretical data by matching one line since none of the reported 
measurements are absolute. These comparisons demonstrate the accuracy of our procedure. In Figs.[4][5]and[6]we also compare 
synthetic absorption spectra of CaH and MgH at the room temperature as well as at T = 1500 K generated using our and 
the UGAMOP line lists. The line intensities in these figures are given by sticks. For a better comparison the UGAMOP line 
intensities are divided by 2 in order to account for the difference in treating the fine (i.e. spin-rotational) splitting. In the case 
of CaH our intensities in the higher frequency range are systematically weaker than that of UGAMOP. Fig. [7] illustrates a 
1500 K absorption spectrum for BeH. 



3.1 Partition function 



We used the calculated energies of all three molecules to generate the partition function values for a range of temperatures 



from to 3000 K. These values were represented by the following polynomial expansion (Vidler & Tennyson (2000l) via a 
least-squares fit: 



log(Z) = ^a < log(T) i 



(8) 



The resulting expansion constants a,i for BeH, 24 MgH, and 40 CaH are listed in Table [H] To a very good approximation the 
partition functions of 24 MgH, 25 MgH, and 26 MgH are in the ratio 2:12:2, which takes account of the increased nuclear spin 
statistical weie ht of 25 MgH. This scaling can be used to obtain Z(T) for the latter two species. 

We have also compared our Z(T) of CaH to the partition function values computed using the UGAMOP energies of CaH. 
Even though the difference in our and UGAMOP energies is large (Table |3|, the partition functions agree within 3% for the 
whole range of T up to 3000 K. 



Irwin ( 1981 1 and Sauval & Tatum ( 1984 1 reported similar expansions of Z(T) for a large set of molecules, which we 
can use to validate our results. For example, for 24 MgH at T = 2000K from Irwin (19811 and Sauval & Tatum (19841 one 
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— 



R(9) 



T 



T 



T 



T 



R(10) 



Figure 3. Theoretical stick absorption spectra of BeH, 24 MgH, and 40 CaH at the room temperature compared to the experimental 
data reported by Shaycst eh et al.| ( |2003^ , [Shaycstch ct al. (2004a|, and |Shayesteh et al.| |2~004b[ l, respectively. The experimental spectra 
are scaled to match the theoretical intensities. 
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Dbserved e 

uculated 

= 2->l 
(9) 

f 

v= 3-^2 
<jf P(4) 


f 

v= l->0 
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- I 
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1 1 r 
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P(19) 



Calculated 



e f 



v 3^2 
P(12) 



v= 4->3 
P(8) 



_L 



1069 1070 1071 1072 1073 1074 

Wavenumber (cm" 1 ) 



1075 



1076 



1077 



© 2012 RAS, MNRAS 000,[T|fT4l 



ExoMol line lists I: The rovibrational spectrum of BeH, MgH and CaH in the X 2 E + state 9 



Table 3. A comparison of the theoretical and experimental spectroscopic constants of "CaH in cm 





UDb. 


Hair- & 


UE)a. — V_j alC . 


11C AMOP 


Ti 


1260.12775 


1260.1275 


0.0003 


1248.1020 


T-2 


2481.99888 


2481.9986 


0.0003 


2459.1800 


T 3 


3665.4141 


3665.4135 


0.0006 


3634.0060 


Ti 


4809.9464 


4809.9453 


0.0011 


4772.1695 


T 5 




5914.8483 




5873.3462 


T 6 




6979.2865 




6936.7620 


T- 




8002.3180 




7957.0594 


T 8 




8982.6564 




8926.4122 


Tg 




9918.3259 




9836.9332 


Tio 




10806.3117 




10680.8451 


Tn 




11642.2418 




11447.2650 


Tl 2 




12420.0739 




12118.8518 


T i3 




13131.7349 




12676.5357 


1X4 




13766.6602 




13121.7989 


Tib 




14311.2476 




13461.1465 


Tie 




14748.4857 




13696.1992 


T 17 




15059.0529 






So 


4.2286902 


4.22868360 


0.66xl0~ 5 




Bj 


4.131722 


4.13171966 


0.23X10" 5 




B 2 


4.0342454 


4.03424554 


-O.OlxlO" 5 




B 3 


3.935887 


3.93589433 


-0.73xl0~ 5 




B 4 


3.836122 


3.83614878 


-2.68X10- 5 




B s 
Be 




3.73452399 
3.63062100 






B- 

B S 




3.52382050 
3.41293766 






Bg 




3.29601845 






Bio 




3.17024576 






B X1 




3.03183969 






B12 




2.87582915 






B13 




2.69557438 






B M 




2.48189083 






B15 




2.22158167 






Bie 




1.89536199 






B17 




1.47691639 







Experimentally derived by Shaycstch ct al. (2004b); Calculated, this work 




4000 6000 8000 10000 12000 
-i 



wavenumbers /cm" 



Figure 4. Absorption spectra of CaH at T=296 K: ExoMol vs UGAMOP 
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2000 4000 6000 8000 10000 12000 

wavenumbers /cm' 1 

Figure 5. Absorption spectra of CaH at T=1500 K: ExoMol vs UGAMOP. 




2000 4000 6000 8000 10000 

wavenumbers /cm" 1 

Figure 6. Absorption spectra of MgH at T=1500 K: ExoMol vs UGAMOP. 




2000 4000 6000 8000 10000 12000 14000 16000 

wavenumbers /cm" 1 



Figure 7. ExoMol absorption spectra of BeH at T=1500 K. 
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Table 4. Extracts from the line list files for MgH, CaH, and BeH. 





1 


2 


3 


4 


5 


6 


7 


8 


9 



91 .048711 


2.142280E-01 


390 031 578 


g 


f 


o 


7 


f 







2.430318E-01 


390 997037 


g 




o 


7 




o 


91 .972295 


3.231341E-01 


4472.543025 


<l 


f 


■-j 


g 


f 


3 


91.993639 


3.614000E-01 


4472.733852 


9 




3 


g 




3 


92.716937 


3.674842E-01 


6900.960666 


10 


f 


5 


9 


£ 


5 


92.734931 


4.063714E-01 


6901.144783 


10 


c 


5 


9 


c 


5 


92.839961 


1.965889E-01 


9918.094815 


14 


f 


8 


13 


f 


8 


92.846033 


2.111897E-01 


9918.247096 


14 


c 


8 


13 


c 


8 


92.884148 


3.860760E-03 


10261. 123179 


1 


c 


1 1 


2 


c 


10 


92.893287 


2.144867E-03 


10261.112833 


1 


f 


1 1 


2 


f 


10 






40 CaH 














1^71 R9fi& 

i-O . j lu^UO 




3673 9fi7'^nn 


9 


f 


3 




f 


3 


1. <J . j U / U O O 


5.424912E-03 


3673 3955X6 


9 




■j 






3 


1 1 1 OQ^Q 
J-U. j. x u ^ 


9 0^71 07P, (V*. 


9400 047709 


2 


f 


2 




£ 


2 


16. 132829 


u. i j u -jyjo hj-\j o 


13147.887247 


3 


f 


13 


2 


f 


13 


16.151200 


5.322792B-03 


2490.108168 


2 


e 


2 


1 


c 


2 


16.154575 


6.453246B-03 


13147.941920 


3 


c 


13 


2 


c 


13 


16.500092 


2.878773B-03 


1268.370045 


2 


f 


1 


1 


f 


1 


16.541842 


5.181792B-03 


1268.432737 


2 


c 


1 


1 


c 


1 


16.887199 


2.782460B-03 


8.434990 


2 


f 





1 


f 





16.930436 


5.008428B-03 


8.499915 


2 


c 





1 


e 









9 BoH 














19.079220 


1.306725B-04 


3896.871469 


1 




2 







2 


19.707403 


9.607023B-05 


1986.417648 


1 




1 







1 


20.327181 


6.337485B-05 


0.000000 


1 















22.576116 


6.097568B-04 


15721.783485 


2 




10 


1 




10 


23.667196 


2.7139S3B-04 


16426.182379 


3 




11 


2 




11 


26.231114 


1.246778B-03 


14731.961064 


2 




9 


1 




9 


28.814926 


1.753481B-03 


13541.079476 


2 




8 


1 




8 


30.852363 


2.056734B-03 


12198.490607 


2 




7 


1 




7 


30.856982 


5.807508B-04 


16449.849576 


4 




11 


3 




11 


32.574171 


2.147658E-03 


10732.929225 


2 




6 


1 




6 



Column Notation 



1 


Pij 


Transition frequency in cm -1 


2 


Aij 


Einstein A coefficient in s — 1 


3 


E" 


lower state energy in cm -1 


4 


N' 


upper state rotational quantum number 


5 


P 


upper state rotationlcss parity 


(i 


v' 


upper state vibrational quantum number 


7 


N" 


lower state rotational quantum number 


8 


p" 


lower state rotationlcss parity 


9 


V 


lower state vibrational quantum number 



Table 5. The completeness and coverage of our line lists for BeH, MgH, and CaH. 





9 BeH 


24 MgH 


25 MgH 


26 MgH 


40 CaH 


I'max (cm 1 ) 


16 400 


10 354 


10 354 


10 355 


15 278 


^max 


11 


11 


11 


11 


19 




59 


60 


60 


61 


74 


Mines 


3968 


6716 


6751 


6754 


26 980 



obtains Z = 659.8 and Z — 783.2, respectively to compare with our value Q/g s ~ 811.0. The factor g s = 2 here is needed to 
account for the nuclear spin statistics absent in the approach used by Irwin ( 1981 \ and Sauval & Tatum ( 1984[ | (see discussion 



Tennyson & Yurchenko (20121). However at room temperature even a small difference in Z(T) becomes important. Our 



value at T = 296 K is 72.8, while Sauval & Tatum ( 19841 give 65.7. These values are also collected in Table [6] together with 
corresponding Z(T) obtained for CaH and BeH. Note a striking difference of about 50% for Z(296K) of CaH between our 
value and that from Sauval & Tatum (1984 1. 
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Table 6. BeH, 24 MgH, and 40 CaH partition function, Z(T), expansion parameter s for Eq. (|8p and sample partition function values. 
These values are compared to Z(T) obtained using parameters from UGAMOP and |Sauval fc Tatum| |l984| l. Please note that our and 
UGAMOP's values are computed using the nuclear statistical factors 8, 2, and 2 for BeH, 24 MgH, and 4u CaH, respectively. 



Parameter 




BeH 


24 MgH 


40 CaH 


Oo 




-3.69(26) 


-4.73(42) 


-2.38(36) 


0,1 




14.63(78) 


16.1(12) 


9.2(11) 


a 2 




-17.95(91) 


-19.8(14) 


-11.5(13) 


«3 




11.64(54) 


13.01(86) 


8.01(75) 


04 




-4.01(18) 


-4.61(28) 


-2.98(25) 


OS 




0.703(30) 


0.836(48) 


0.564(42) 


(-Hi 




-0.0490(21) 


-0.0606(33) 


-0.0423(29) 






BeH 


24 MgH 


40 CaH 


9ns 




8 


2 


2 


Z(296K) x 9ns 


SfcT" 


93.9 


129.5 


93.9 


Z(296K) 


UGAMOP 




147.8 


195.9 


Z(296K) 


This work 


335.0 


145.6 


199.9 


Z(2000K) xg ns 


S&T a 


2276.9 


1566.4 


2276.9 


Z(2000K) 


UGAMOP 




6575.5 


2451.7 


Z(2000K) 


This work 


3019.5 


1622.4 


2253.5 



a |Sauval fc Tatum| |l984] 
4 CONCLUSIONS 

Our potentials are highly accurate for low vibrational levels but there is undoubtedly some loss of accuracy for higher 
vibrational states given the lack of experimental data to constrain the curves. For BeH, high accuracy should be achieved for 
all transitions involving states up to v = 11. For MgH, we can be confident of sub-wavenumber accuracy for v > 4 levels as 
these higher levels can be characterised using electronic transitions. For CaH, it is not possible to test any of the vibrational 
states higher than v = 4 since there is no experimental data. 

Line lists for the rotation-vibration transitions within the ground states of the BeH, 24 MgH, 25 MgH and 26 MgH, and 
40 CaH molecules are presented which should be accurate over an extended temperature range. These line lists were computed 
using potential energy curves fitted to experimental data and ab initio dipole moment curves obtained at high levels of theory. 
A complete line list for each of these molecules can be downloaded from the CDS, via ftp : //cdsarc .u-strasbg. f r /pub/cats/ 
J/MNRAS/ , or via http : //cdsarc .u-strasbg. fr/viz-bin/qcat?J/MNRAS/ The line lists together with auxiliary data including 
the potential parameters, dipole moment functions, and theoretical energy levels can be also obtained at www.exomol.com. 
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